Pressure-Driven Metal-Insulator Transition in Hematite from Dynamical Mean-Field 
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The Local Density Approximation combined with Dynamical Mean-Field Theory (LDA+DMFT 
method) is applied to the study of the paramagnetic and magnetically ordered phases of hematite 
Fe2 03 as a function of volume. As the volume is decreased, a simultaneous 1st order insulator- 
metal and high-spin to low-spin transition occurs close to the experimental value of the critical 
volume. The high-spin insulating phase is destroyed by a progressive reduction of the charge gap 
with increasing pressure, upon closing of which the high spin phase becomes unstable. We conclude 
that the transition in Fe2 03 at «50 GPa can be described as an electronically driven volume collapse. 

PACS numbers: 74.25.Jb, 71.27.+a, 71.30.+h, 91.60.Gf 



Metal-insulator transition (MIT) is one of the central 
topics in the physics of strongly correlated electron ma- 
terials. Despite considerable progress in the theoretical 
understanding of this phenomenon in the past decade, 
we are only starting to uncover the multitude of possi- 
ble transition scenarios in real materials. An important 
instance of the MIT is the pressure driven transition ac- 
companied by a change of the local spin state (high spin 
(HS) to low spin (LS) transition) seen in MnO [jj, BiFe0 3 
[2| or Fe 2 03. The understanding of the pressure-driven 
HS-LS transition and its relationship to eventual MIT 
and structural and/or volume changes is relevant to a 
broader class of oxides - often with important geological 
implications - such as FeO, CoO or Fe304- 

In this Letter we study the spin and metal-insulator 
transition in hematite (a-Fe203) under pressure using 
the numerical LDA+DMFT approach [IJ. At ambient 
conditions, hematite is an antiferromagnctic insulator 
(T/v=956 K) with the corundum structure The iron 
ions having a formal Fe 3+ valency with five d electrons 
give rise to a local HS state. Photoemission studies [!, 0] 
classified hematite as charge-transfer insulator. A gap of 
2.0-2.7 eV was inferred from the electrical conductivity 
data 0- Under pressure, a first-order phase transition is 
observed at approximately 50 GPa (82% of the equilib- 
rium volume) at which the specific volume decreases by 
almost 10% and the crystal-lattice symmetry is reduced 
(a structure of the RI12O3-II type is formed) [1, [tj [Tol ] . 
The high-pressure phase is characterized by a metallic 
conductivity and absence of both the long-range mag- 
netic order and the HS local moment ||. Badro et al. 
showed that the structural transition actually precedes 
the electronic transition, which is nevertheless accompa- 
nied by a sizable reduction of the bond lengths [11 1 . 

The investigations of the electronic structure of 
hematite using various theoretical approaches produced 
different (and frequently controversial) results. Simi- 
lar to many transition metal oxides, calculations within 



the local spin density (LSDA) and generalized gradient 
(GGA) approximations failed to reproduce the ground 
state characteristics under ambient conditions, underes- 
timating values of the magnetic moment of iron and of 
the energy gap [12j . Within GGA, the structural transi- 
tion together with a transition into the metallic ferromag- 
netic state with Fe ions in the LS state was predicted at a 
pressure of 14 GPa [131 ]. which is about one fourth of the 
experimental value. The studies using the LSDA+U ap- 
proach, which treats the local Coulomb interaction on the 
static mean-field level 13j lead to improved values of the 
magnetic moment and the gap in the HS phase, but failed 
to capture the metal-insulator transition. Like for some 
other oxides, LSDA+U improves the description of the 
HS low-pressure phase, while LSDA appears more appro- 
priate for the description of the LS high-pressure phase, 
none of them covering the full pressure range. Thus, the 
DMFT method, which can capture the evolution from 
strong to weak local correlations is a promising tool for 
investigating the physics in the transition regime, as was 
shown for MnO Q. The only available DMFT study 
of Fe203 to-date was limited to the t2 g band manifold 
finding a low-spin metallic phase at high pressure [1 51 ] - 

In this work we use a Hamiltonian including explicitly 
all O-p and Fe-d orbitals, which allows us to address the 
following questions, i) Is there an electronic transition 
without invoking a change in crystal structure? ii) Is 
there a simultaneous metal-insulator and local moment 
transition as found in the isoelectronic MnO |l4|? iii) 
What is the nature of the transition and how is it affected 
by temperature and long-range magnetic order? We use 
an implementation of LDA+DMFT, previously applied 
to NiO and described in Rcf. which combines the 
first-principles multi-band Hubbard Hamiltonian 

H = £)&tH fc fi* _ £dc {n d )N d + £ h d V dd h d 

k ieFe 

with DMFT equations. Here, H fe is the LDA Hamil- 
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FIG. 1: Orbital occupancies (left panel) at various specific 
volumes in PM state at T=580 K. The thick line marks the 
average d occupancy per orbital. Corresponding values of the 
screened (diamonds) and instantaneous (circle) local moment 
are shown in the right panel. 



tonian matrix (38x38 plus spin degeneracy; two Fe2 03 
units) constructed on the basis of Fe-d and O-p Wan- 
nier orbitals represented numerically on a uniform k- 
mesh in the first Brillouin zone [T3|. The second term is 
the double-counting correction amounting to a constant 
shift applied to Fe-d site energies. The last term is the 
two-particle interaction at the Fe sites in the density- 
density approximation, parameterized with [7=6.8 eV 
and J=0.86 eV [T3|. (Fitting the experimental photoe- 
mission spectra yields U of 8 eV according to Ref. @ or 
7 eV according to Ref. For the double-counting 

correction Edc we have adopted the prescription edc = 
(N OI b — l){7n<2 of Ref. Il6|, but taking the self-consistent 
DMFT value for the average occupancy fid rather than 
its LDA value. This approach, which can be viewed as a 
poor man's charge self-consistency loop, is necessary for 
stabilization of the ambient pressure phase as discussed 
later. For the solution of the <i-shell impurity problem we 
have employed a continuous time quantum Mote-Carlo 
(CT-QMC) solver [joj . The calculations were performed 
in the corundum a-Fe203 structure with rhombohedral 
(space-group R3c) unit cell containing two formula units. 
Throughout the paper the compressed phase is character- 
ized by the volume relative to the experimental ambient- 
pressure structure with lattice parameter ao=5.361 A, 
rhombohedral angle a=55°17', and atomic coordinates 
Fe(0,0,0.355) and 0(0.300,0,0.25). Volume reduction is 
achieved by scaling the lattice constant. 

In Fig. Q] we show the evolution of the orbital occu- 
pancies and Fe local moment under compression in the 
absence of long-range magnetic order. The local moment 
is measured as the square root of i) {ml) (instantaneous 
moment) or ii) ^ G?r(m z (r)r7i z (0)} (screened moment), 
where (3 is the inverse temperature. At the ambient con- 
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FIG. 2: (Color online) Evolution of the single-particle spec- 
tra with pressure obtain in the PM state at T=580 K. The 
HS solutions are shown in the left and LS solutions in the 
right panel. All curves are normalized to one. The left 
inset presents comparison of experimental PES 6] and IPS 
[22| data (points) to the calculated Lorentzian broadened d- 
spectral density (solid line). The right inset shows a compar- 
ison of the v=0.8 DMFT (smooth lines) and non-interacting 
LDA spectra (jagged lines): total d spectral density (dark; 
black), total O-p spectral density (light; blue). 



ditions (volume Vq) we find a HS solution with both def- 
initions of local moment yielding essentially the value of 
4.6 fiB- Upon compression the HS phase survives down 
to 0.8Vb- Below this volume only a LS phase with sub- 
stantially reduced local moment is stable. The difference 
between its screened and instantaneous values indicates 
that the LS moment cannot be viewed a rigid object. The 
HS and LS solutions coexist in the range 0.8 — 0.85Vo- 
The Fe-d orbital occupancies fit into the ionic picture of 
a fully spin-polarized d-shell in the HS state, and an or- 
bitally polarized d-shell in the LS state. Nevertheless, 
the substantial e£ population in the LS phase reflects a 
sizable covalent p— d bonding (hybridization between O-p 
and Fe-e ff orbitals). 

In Fig. [2] we show the corresponding single-particle 
spectral densities. Note the difference of the O-p spec- 
tral weight above the chemical potential between HS and 
LS solutions. In the LS state the e^-derived bands are 
above the chemical potential, carrying some admixture 
of O-p character with them, which leads a charge trans- 
fer from the O-p orbitals and an increase of fid (see Fig. 
[T]). Relaxation of the electronic states, which compen- 
sates this charge redistribution in real material, is mim- 
icked by the self-consistency of the double counting term 
in the present work. Comparison to the photoemission 
and inverse photoemission data (see left inset of Fig. ^ 
reveals a fair agreement between the theory and experi- 
ment. The computed ambient pressure HS phase is found 
to be an insulator with a charge gap of 2.6 eV. Upon com- 
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FIG. 3: Spin-polarized spectral Fe-d spectral density at ambi- 
ent pressure for 580 K, 1160 K and 1450 K from the top to the 
bottom. The left inset compares the same densities averaged 
over spin. The right inset shows the staggered magnetiza- 
tion vs temperature curves for v=1.0 (circles; black), v=0.9 
(squares; red), v=0.85 (diamonds; green), and v=0.8 (trian- 
gles; blue). At empty symbols only the LS solution exists. 
The lines depict mean-field fits. 



pression the gap gradually shrinks due to broadening of 
the bands and increasing of crystal-field splitting. At the 
high-pressure end of the HS stability region (O.8V0) the 
gap is reduced to a fraction of an eV. The spectral density 
in the LS phase closely resembles the non-interacting one 
(see right inset of Fig. Unlike in the HS phase, there 
is some orbital polarization within the t2 g {— aig+el) or- 
bitals, split due to distortion of FeC"6 octahedra. 

Next, we address the role of magnetic long-range order. 
Unlike LDA or LDA+U, for which static magnetic order 
is essential to capture a HS insulator, in the present ap- 
proach magnetic ordering can take place independently 
from the gap and local moment formation, as observed in 
many materials. Using a symmetry that allows for the ex- 
perimental AFM order we observe the following behavior. 
At room temperature the HS solutions are magnetically 
ordered, while all LS solution remain paramagnetic. In 
Fig. [3]we show the spin-polarized Fe-ci spectral densities 
at several temperatures. Clearly, the charge gap and the 
main spectral features are insensitive to the presence of 
AFM order. The magnetization curves at reduced vol- 
ume reflect the increase of the exchange coupling with 
shortened atomic distances. 

Our observations are summarized in the phase diagram 
shown in Fig. |U We distinguish three phases: i) para- 
magnetic LS, ii) antifcrromagnetic HS and iii) param- 
agnetic HS, separated by the line of Neel temperature 
Tjv and the coexistence region of the HS and LS phases. 
The evaluation of TV is affected by two approximations, 
both favoring the ordered phase. First, the lack of non- 
local fluctuations inherent to DMFT and second, the ap- 
proximate form of the local Coulomb interaction, limited 
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FIG. 4: Calculated V-T phase diagram of hematite. The 
symbols mark the actual computation results: HS state (cir- 
cle), LS state (diamond). Empty diamonds mark LS solu- 
tions stable only if PM symmetry is constrained, dotted cir- 
cles mark HS solutions stable only if AFM order is allowed. 
The black line shows the Neel temperature estimated from 
fits to v=1.0, 0.9 and 0.85 data. The shaded area marks the 
estimated HS/LS coexistence region with PM constraint, the 
dashed lines mark the coexistence regime when AFM order is 
allowed. 



to Ising terms, thus rendering the HS local state two- 
(S z = ±5/2) instead of sixfold degenerate (S = 5/2). 
The latter approximation is also reflected in the shape of 
the magnetization curves, which can be well fitted with 
a mean- field theory based on the Brillouin function B 1 / 2 
(see the inset of Fig. [3]). Given these approximations 
the calculated TV of 1600 K appears to be in reasonable 
agreement with the experimental value of 956 K. 

The effect of AFM order on the HS-LS transition is to 
slightly shift the boundaries of the coexistence region in 
favor of the HS phase. In particular, at the lowest stud- 
ied temperatures we find the LS solution at 0.85Vo to be 
unstable when AFM order is allowed. Also we find the 
HS solutions at (0.8F , H60 K) and (0.825V , 1660 K) 
stabilized by AFM order. Interestingly, in the vicinity 
of the HS /LS boundary we find the magnetic order with 
increasing temperature not to be destroyed by magnetic 
fluctuations, but via a collapse of the constituting mo- 
ments. 

We propose two scenarios how the HS state can be 
destabilized and discuss them in the context of the 



present results, the MnO results of Ref.ll4l and the model 
calculations of Ref. [2l|. We start with a simple estimate 



of the charge gap considering the ionic energies 

E e w E{d 6 )+E(d 4 )~2E(d 5 )-W = U+4J-2A-W, (1) 

where 2A stands for the difference of ^2g(= a ig+ e g) an d 
energies (as in Ref. ,21), W for the average of the 
e g and t 2g bandwidths, E(d 5 ) is the energy of the HS 
ionic groundstate while E(d 6 ) and E(d 4 ) the ionic en- 
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ergies of the state with one more and one less electron, 
respectively. Similarly one can compare the energies of 
the ionic HS and LS states to get 



E{LS) - E(HS) = 10 J - 4A. 



(2) 



The HS state becomes unstable when either ([T]) or ([2]) 
becomes negative. Thus depending on the parameters 
the HS state can be destroyed due to gap closing or due 
to a local state transition. In charge-transfer compounds 
such as hematite E g is substantially reduced by the pres- 
ence of the ligand bands. We found that a similar simple 
estimate for a 2-band model provides a fairly accurate 
description of the phase boundaries found in Ref. |2U 

Comparing the present results with the earlier study 
of isoelectronic MnO 14[ there is an obvious similarity 
in that the HS phase is an insulator while the LS phase 
is metallic. However, there are two important differences 
in how the transition happens. First, in MnO the tran- 
sition is found to proceed through a continuous sequence 
of intermediate states while in hematite the transition is 
first order with a coexistence region of HS and LS so- 
lutions. Second, we find the gap to almost close at the 
high- pressure end of the HS stability region in Fe2C>3 
while in MnO the transition starts with the appearance 
of an in-gap spectral density. The latter observation leads 
us to suggest that the HS-LS/insulator-metal transition 
in hematite proceeds through the gap closing mechanism 
while in MnO a local state transition takes place as was 
argued in Ref. 1 141 The gap closing picture of the tran- 



sition in hematite is also supported by the shape of the 
HS/LS phase boundary, showing that HS state can be 
destroyed by thermal fluctuations, contrary to the local 
state transition picture where the HS state is expected to 
be favored at high temperature due to its higher entropy. 

Finally, we briefly address the reported structure and 
volume changes 0,0, [Hj]. The corundum and Rh 2 03-II 
structures have similar local Fe environment and transi- 
tions between the two were reported in other materials 
without significant changes of electronic structure. In 
this study we do not attempt to distinguish between 
the two structural types. Finding the electronic HS- 
LS/insulator-metal transition in the high- volume phase 
shows that the electronic transition is not implicated by 
the structural transition or volume collapse. Conversely 
the electronic transition provides a natural explanation 
of the volume collapse as a consequence of emptying of 
the anti-bonding e a g bands in the LS phase and corre- 
sponding strengthening of the Fe-0 bonds. 

In summary, we have observed a simultaneous HS- 
LS/insulator-metal transition in the corundum phase of 
hematite close to the experimental volume in the numer- 
ical results obtained by LDA+DMFT. The long-range 
magnetic order is found to play a minor role in the tran- 
sition. We have discussed two possible scenarios, the 
gap closing and the local state transition, and concluded 
that the transition in hematite, unlike in MnO, proceeds 



through the gap closing mechanism. Our results are con- 
sistent with the picture of an electronically driven volume 
collapse. 
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